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Abstract: We study the role of scale dependence in the Buchert averaging method, using 
the flat Lemaitre-Tolman-Bondi model as a testing ground. Within this model, a single 
averaging scale gives too coarse predictions, but by replacing it with the distance of the 
objects R{z) for each redshift z, we find an precision at z < 2 in the averaged 

luminosity and angular diameter distances compared to their exact expressions. At low 
redshifts, we show the improvement for generic inhomogeneity profiles, and our numerical 
computations further verify it up to redshifts 2; ~ 2. At higher redshifts, the method breaks 
down due to its inability to capture the time evolution of the inhomogeneities. We also 
demonstrate that the running smoothing scale R{z) can mimic acceleration, suggesting 
it could be at least as important as the backreaction in explaining dark energy as an 
inhomogeneity induced illusion. 



Keywords: Dark Energy, Inhomogeneous Cosmological Models, Cosmology, Gravitation 



* E-mail: teppo.mattsson@helsinki.fi 
t E-mail: maria.ronkainen@helsinki.fi 



Contents 



i 



Introduction 

The Buchert equations for the flat LTB model 

|2.1| Observations in the LTB model 

|2.2| Observations in the averaged LTB model 

Observables in scale dependent Buchert averaging 

H Four levels of coarse graining 

Analytic considerations — small z behavior 
|3.3| Acceleration without backreaction 



s 



3.5 



3.4 Periodic inhomogeneities as a toy model for structure 
Discussion of the results 

R[z) = R — Single scale 

R{z) = f{z) — Running scale with averaged geodesies 
3.5.3 Riz) = rLTB(-z) — Running scale with exact geodesies 



3.5.1 



3.5.2 



10 
11 
14 
17 
17 
17 
19 



El Conclusions 



19 



1. Introduction 

The current cosmological observations seem to get the simplest and rather concordant 
interpretation in the homogeneous and isotropic expanding universe models, with late- 
time acceleration starting around the redshift z ~ 0.6 [1-3]. The acceleration is usually 
seen as an evidence for dark energy, most often in the form of a cosmological constant or 
vacuum energy. However, the enormous fine-tuning needed to explain both the size and 
the timing of such an energy component has raised serious doubts about its correctness 
and thus justified the search for alternatives [4-9]. 

Perhaps the most natural alternative explanation so far has arisen from the inhomo- 
geneous cosmological models. The point in these models is that suitable inhomogeneities 
can have a similar effect on the observations of light as accelerating expansion in the ho- 
mogeneous models [10-13]. Although gained more popularity only recently, the actual idea 
is not a new one (see [14, 15]). Indeed, already the pioneers of cosmology were careful to 
point out the potential inadequacy of the simplest homogeneous models in describing the 
real universe [16-18]. 

Two conceptually rather different kind of inhomogeneities have been proposed as the 
culprit for the apparent acceleration. Firstly, the non-perturbative effects of the well- 
established lumpiness of galaxies and galaxy clusters are still unknown, and could poten- 
tially mimic acceleration [19]. A virtue in this scenario is that it would connect the growth 
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of nonlinear structure with the start of the acceleration era and thus solve the coincidence 
problem [20]. Secondly, there are inhomogeneities also on scales larger than the observed 
clustering. Indeed, the increasing accuracy of the cosmological observations has revealed 
larger and larger voids [21, 22]. For these kind of smooth large scale inhomogeneities, there 
are exactly solvable models, such as the LTB class of solutions, which can mimic acceler- 
ation without dark energy [11, 13]. However, due to the complexity of the Einstein field 
equations, there are no exact solutions for the small scale lumpiness of the universe. Thus, 
some level of coarse graining has to be introduced, which has usually been done in the 
form of an averaging method [23]. One of the most popular methods in cosmology is the 
Buchert averaging [24, 25], which is also the one considered in this work. 

The conventional way to apply the Buchert formalism is to average over a single do- 
main, larger than the supposed scale of statistical homogeneity in the galaxy distribution. 
Consequently, it has been speculated that the averaging method fails if there are inho- 
mogeneities at large scales as well [19]. Motivated by the recently observed large voids 
and superclusters, we try to extend the averaging method to work with large scale inho- 
mogeneities by going beyond the single scale approximation. In fact, the point we want 
to bring out is that it would be physically reasonable to replace the single scale R by the 
distance of the objects R{z) for each redshift z, since the distance the observed light travels 
depends on how far the object is. 

A commonly presented conjecture in cosmological averaging is that the backreaction 
of inhomogeneities causes the acceleration of the average expansion and could thus account 
for the observations [19, 25, 26]. However, as demonstrated in this work, by promoting the 
single averaging scale R to the redshift dependent function R{z), it is possible to mimic 
acceleration even in the absence of backreaction. Indeed, at least within the flat LTB 
model, this extension gives rise to a definite improvement in accuracy of the averaged 
luminosity and angular diameter distances compared to the single scale case, suggesting 
that the averaging method can also be utilized for large scale inhomogeneities. Naturally, 
in the case a single scale would suffice, the running scale approach reduces to give the 
same predictions as the conventional single scale approach. An additional virtue is the 
computational simplicity of this generalization. 

The paper is organized as follows. In Sect. ^, we derive the observable distance-redshift 
relations for the averaged LTB model with flat spatial sections. The principal results are 
presented in Sect. ^, where we compare these relations with their exact counterparts, using 
various implementations of the running averaging scale R{z). For general inhomogeneity 
profiles, we utilize power series to make an analytic comparison at low redshifts. At higher 
redshifts, we employ numerical computations for two explicit profiles: a bubble inhomo- 
geneity that fits the supernova observations and periodic inhomogeneities as a toy model 
for structure. The results of the different cases are discussed in Sect. 3^. Finally, Sect. § 
contains our conclusions. 

2. The Buchert equations for the flat LTB model 

In this section, we calculate the Buchert equations for the spatially flat, spherically sym- 
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metric LTB model with pressureless matter as the only source. The backreaction of the 
model vanishes identically [26] and hence the only difference from the homogeneous and 
flat matter dominated FRW case is the scale dependence of the averaged quantities. This 
property makes the model especially useful in extracting the effect of the averaging scale on 
the observable quantities, such as the relation of redshift to the angular diameter distance 
and to the luminosity distance. 

2.1 Observations in the LTB model 

The line-element of the spatially flat LTB model with the spatial origin at the symmetry 
center reads as 

ds^ = -dt^ + {A'{r, t)fdr'^ + A^{r, t) {d0^ + sin^ edip^) , (2.1) 

where A{r, t) is the scale function having both temporal and spatial dependence, and we 
use the following shorthand notations for the partial derivatives: ' = ^ and ' = ^. This 
metric was first studied by Lemaitre [16], Tolman [17] and Bondi [18]; later, it has been 
used in various astronomical and cosmological contexts [14, 15]. Although commonly called 
a toy model, the metric ( |2.1| ) is an exact solution of the Einstein equations and the perfectly 
homogeneous FRW model is only a special case of it, obtained in the limit: A(r, t) — > a{t)r, 
where a{t) is the FRW scale factor. Our notation and parametrization follows Ref. [27]. 

The Einstein equations for the metric ( |2.lD reduce to the generalized Friedmann equa- 
tion 

where H{r,t) = A{r,t) / A{r,t) is the LTB version of the Hubble function, AQ{r) = j4(r, tg) 
is the scale function at a reference time to) HQ{r) = H{r,to) is the position dependent 
Hubble constant, and to the time evolution equation of the matter density 

Al{r)A',{r) 



1A^{r)%{r) 
^ 3A',{r)Ho{r) 



A'^{r,t)A'{r,t) 



The integration of Eq. (|2.2D w.r.t. time yields 



A{r,t) 



2/3 



Ao{r) 

Substituting Eq. ( |2.4D in Eq. ( |2.2| ) gives the time evolution of the Hubble function: 

^o(r-) 



(2.3) 



(2.4) 



H{r,t) = . (2.5) 

l + ^(t-to) 

The choice of AQ{r) represents a coordinate freedom, similar to the normalization of a(to) 
in the FRW case; here we set ^o(^) = Consequently, the flat LTB model is uniquely 
determined by the free function Hq^v). Plugging Eq. ( p.4| ) in Eq. ( |2.3| ) then gives the 
explicit time dependence of the matter distribution: 

' 87rG[l + ^(t-to)][r(i-to)W + (l + ^(*-*o))] ' 
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Eqs. (2^) and show that the inhomogeneities in the flat LTB model correspond 
to "decaying modes" [28], i.e. the model evolves towards homogeneity at late times. This 
is just the opposite to the expected structure formation in the real universe, but for low- 
redshift observations, this should not make a difference. In the nonflat LTB model with 
simultaneous Big Bang, the inhomogeneities correspond to "growing modes" and in this 
respect, that model would be more realistic. However, as the backreaction does not vanish 
in the nonflat model, it would be harder to distinguish the effect of the running smoothing 



scale from the backreaction effects, which is why we use the flat LTB metric (2T). The 
price to pay is that, for most of the inhomogeneity profiles HQ{r), we have to restrict our 
considerations to observations at low redshifts (z < 2). 

To study observable properties of light, we need the pair of radial geodesic equations 
determining the relations between the coordinates and the observable redshift, t{z) and 
r{z), given by [18] 

dt _ -A'{r,t) 
dz 



dr 
dz 

To calculate the right hand sides of Eqs. (|2^ 
of the scale function (p.4|): 



{l + z)A'{r, t) 
1 

(l + z)i'(r, t) 



(2.7) 



(2.8) 



and (|2.8|), we need the following derivatives 



A'ir,t) 



l + -(t 
2^ 



to)^o(r) + r{t - to)H^{r) 



l + -(t 
2^ 



to)Ho{r) 



(2.9) 



i'(r,t) 



to)Ho{r) 



Hoir) 



to)Ho{r) + 



+rH^ir)(^l + {t-to)Ho{r)^ 



(2.10) 



The relation of the redshift to the energy flux F, or the luminosity-distance, defined 
as di = y^L/AnF with L the total power radiated by the source, as well as to the angular- 
diameter distance dA are given by [29] 



dUz) = {l + zfAir{z),tiz)) 



(2.11) 



dA{z) = A{r{z),t{z)) 



(2.12) 



As the relations t{z) and r(z) are determined by Eqs. ( |2.7D , (|2.8|), (^), ( |2.10| ), and the scale 
function A{r, t) by Eq. ( |2.4| ), using Eqs. (2.11) and ( |2.12| ), one can compute the observables 
diiz) and dA{z) for a given z. In Sect. ^, we compare the exact forms of these observable 
relations to their Buchert-averaged counterparts with various implementations of the scale 
dependence and different inhomogeneity profiles HQ{r). For this, we need to first apply the 



Buchert formalism to the flat LTB metric (2.1 



- 4 - 



2.2 Observations in the averaged LTB model 

To construct a coarse grained description of the inhomogeneous universe, one usually has 
to average dynamical quantities of Einstein's gravitation theory (see [23] for a review). It 
seems physically more correct to first calculate the Einstein field G(g) for the exact metric 
g and only then average (G(g)), than to calculate the Einstein field for the averaged metric 
G((g)). The reason is that the Einstein field is more closely related to physical quantities 
whereas the metric corresponds to gravitational potentials, whose derivatives determine 
the physics. Since in Einstein's gravity the field G depends nonlinearly on the metric g, 
its evaluation does not commute with averaging: (G(g)) / G((g)). Hence the issue is not 
only a conceptual one, but in general leads to physically different predictions. However, 
in the absence of nonlinear inhomogeneities the two approaches lead to identical results. 
In fact, the standard model of cosmology builds on the assumption that the undoubtedly 
existent, intense small scale lumpiness has no cosmological significance. In any case, the 
inadequacy of the standard model to explain the observations without a severely fine-tuned 
cosmological constant should, at the very least, justify the more thorough considerations 
of this assumption. 

By averaging the scalar part of the Einstein equations in the above explained order, one 
arrives at the Buchert equations describing the averaged dynamics of a general irrotational 
dust universe [24]: 

3^ = -A'kG{p)v + Qv (2.13) 

2,1 =SitG{p)v-\{^''^R)v-\Qv (2.14) 
\av J 2 2 

^Ap)v = -S^{p)v , (2.15) 
at ax> 

where p is the matter density and the difference between the Buchert acceleration equation 



( 2.13 ) and its homogeneous FRW counterpart is known as the backreaction 

Qvit) ^ l{{e^)v - {0)1) - {a'^''c7,,)v , (2.16) 
where (J^^Gnu > represents the shear, 



av{t) ^ ( Vdet[g,,]d^x \ ~' 

V Iv \/det[5jj(t = to)]d'^x J 

is the averaged scale factor, D is the averaging domain, gij is the spatial part of the metric, 
^^^R is the curvature scalar of the t = const, spatial hypersurfaces, 9 = V^u^ is the 
expansion scalar and the spatial average of a scalar S is defined as 



{S)v{t) = m.^-ljy=^ . (2.18) 

^/(let[gij]d^x 



The average expansion accelerates if the right hand side of Eq. ( 2.13 ) is positive; 



this can be achieved by having large enough variance of the expansion rate although it is 
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partially counterbalanced by the average shear. The variance gets large when contracting 
{6 < 0) and expanding (9 > 0) regions coexist, and in fact the average acceleration has 
been sometimes connected to gravitational collapse [19,30,31]. Anyhow, as demonstrated 
in Ref. [26], a globally expanding dust universe can have average acceleration as well. 

In the context of the averaged universe, whose dynamics is given by Eqs. ( |2.13| ), ( |2.1^ ) 
and ( p.l5| ), it is natural to further assume that the average metric takes the FRW form: 



-dr + a^{t) 



1 — k-Dr^ 



(2.19) 



as was recently highlighted by Paranjape and Singh [32] (see also [25]). Although the 
form of the metric ( 2.1S| ) is the same as in the perfectly homogeneous universe, the time 
evolution of the scale factor axi{t) and the spatial curvature kx> are in general different 
from any FRW model. 

Recalling the definition of the shear tensor 



7t(V^Ui. + V ) 



(2.20) 



where is the four-velocity of the dust, we obtain the shear and the expansion scalars of 
the LTB metric in the coordinates of Eq. ( |2.1| ): 

2 



A{r,t) 




A'{r,t) 
A'{r, t) 



2H{r,t) + Hr{r,t) 



(2.21) 



(2.22) 



A{r,t) A'{r,t) 

Using the field equation ( |2.2D , one finds that the two Hubble functions H{r, t) and Hr{r, t), 
defined via Eq. ( 2.22| ), are related at the reference time t = to as 

Hr{r,to) = Ho{r)+rH'o{r) , (2.23) 

so that the expressions for the shear and the expansion rate simplify to: 



to) = -{rH',{r)f 



(2.24) 
(2.25) 

(2.26) 

Now we have all the ingredients to apply the Buchert equations ( 2.13| ), ( p. 14 ) and ( |2.15| ) 
to the flat LTB metric (2.1). For symmetry and simplicity, we only consider averages over 
a spherical domain, denoted TZ, of radius R centered at the origin. Hence, plugging Eqs. 
( 2.24 ), ( |2.25| ) and ( |2.26| ) into the expression of the backreaction ( p.l6| ) yields 



e{t = to) = 3Ho{r) + rH'Q{r) . 
The integration measure for the metric ( |2.1| ) at the t = to hypersurface reads as 

^det[gij{t = to)]d'^x = sinOdrdOdip . 



_3_ 



Qnito) = - 



{rH'Q{r)f)r'^dr 



3 

{{3Ho{r) + rHl,{r)fr^dr + 



_3_ 



R 



{3Ho{r) +rH[^{r)ydr 



(2.27) 
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where we have also used the spatial average as defined in Eq. ( 2.18| ). Rewriting the terms 
of the integrands in Eq. (2.27) as total derivatives, gives 



Qnito 



3 

i?3 



R 



d_ 

dr 



3 



R 



d 
dr 



{r^Ho{r))dr 



(2.28) 



Since the reference time t = to is completely general, we can conclude on the grounds 
of Eq. (|]2|), that 

Qn{t) = (2.29) 

for any time coordinate t and averaging radius R. Altogether, we have verified that the 
backreaction Qn vanishes identically for the flat matter dominated LTB model, as was 
already shown in [26]. As an aside, it does not vanish for an integration domain of arbitrary 
shape; we have checked numerically that for various profiles of HQ^r), the backreaction is 
nonzero e.g. over a cubic region. 

Moreover, transforming the coordinates of the t = to hypersurface shows that the 
spatial metric 



da^ 



dr 



dr^ + Al{r){de'^ + sin^ Odip"^) = dAj + Al^dO^ + sin^ edip"^) 



(2.30) 



reduces to the flat Euclidean form, for which the Ricci scalar vanishes and hence we have 
((3)i?)p = in Eq. ( |2l4| ). 

Since Qn = and ((3)/?)^ = 0, the Buchert equations plsD , ( ^ ) and (^) reduce 
to the corresponding FRW equations for the averaged scale factor an{t) and the averaged 
matter density {p)ti- 

-{p)n , (2.31) 



an 
an 



an 

2 



3 

8^G 



{p)n , 



Tr.{p)n = -3—{p)n ■ 
at a-ji 



(2.32) 
(2.33) 



The only difference in the Eqs. ( |2.31| ), ( |2.32| ) and ( |2.33| ) compared to the FRW equations is 
the dependence on the averaging scale. Therefore, they also have the Friedmann solution, 
an{t) = {.i/to{R))'^^^ with tQ{R) as the scale dependent age of the universe, and the template 
metric ( ^.191) reduces to: 

ds^ = -dt^ + {t/to{R)f^ [dr^ + r'^{de'^ + sin^ Od^'^)] . (2.34) 

One can now apply the metric (|2.34| ) and the averaged equations ( ^.3lD , ( ^.32|) , (|2.33| ) 
to calculate the distance-redshift relations (|2.11| ) and ( 2.12 ). By evaluating the expectation 
value of the expansion scalar, 

jT^dt^/det[gij\d!^x _ dv 



(0) 



V 



Iv \/det\g^d^3 



3— 



(2.35) 
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we obtain a simple relation between the averaged Hubble constant TLniio) = O7^(^o)/o7^(^o) 
and the LTB Hubble function (12? 



nn{to) = \{eo)n = (Hoir) + lrH'oir))n = ^{r^Ho{r))dr = Fo(i?) , (2.36) 



where we have also used Eq. ( p. 251) . Using A{r,t) = afi{t)r with Eqs. (2^), (P^), (|2.3 
and (|]3|), we obtain the averaged distance-redshift relation 

1 

where the averaged Hubble constant TCizito) has been eliminated with the help of Eq. 
( 12:361) . Utilizing Eqs. pll| ), ( ^ ) and { ^l7\ ), the observables can be written in the final 



fiz) = 2H,\R){l-^f=] , (2.37) 



form: 



dL{z) = (1 + zfan{t{z))f{z) = 2H^\R{z)){l + z)[l- ) , (2.38) 



dA{z) = an{t{z))f{z) = '^^l^f^^^^ (l - ^jL=) , (2.39) 

where instead of a single scale i?, we have allowed for a running averaging scale R{z). 
Although then the averaged equations ( p.31| ), ( p.32| ) and ( |2.33| ) are not anymore satisfied 
identically, the equations hold true for each scale R{z) separately. In practice, this means 
that we are considering a different FRW model for each redshift. 

The consequences of the redshift dependent averaging scale in the observable relations 



( 2.38 ) and (|2.3g| ) constitute the subject of the following section. 



3. Observables in scale dependent Buchert averaging 

The conventional way to apply the Buchert equations is to choose a single averaging domain 
P, appropriate for the physical system in question. This approximation has been justified 
qualitatively on the basis of the observed statistical homogeneity and isotropy [19]: For 
example, when taking the averages over an observer-centered ball of radius R, its size would 
have to be large enough compared to the small scale lumpiness. As long as this condition 
is fulfilled — the argumentation goes — the actual value of R would be irrelevant. 

However, when applied to large scale inhomogeneities, the validity to use only a single 
averaging domain has not been examined carefully before, though there have been some 
hints towards the idea of multiple averaging scales [25,32]. In the following, we consider 
this quantitatively with the explicit toy model of Sect. ^, for which it was shown that after 
averaging, the scale dependence remains the only difference from the homogeneous and flat 
matter dominated FRW case. 

To make the treatment physically reasonable, we calculate observables, such as the 
luminosity-distance ^^(z) and the angular-diameter distance dA{z)-, given by the appropri- 
ate averaged expressions ( |2.38| ) and ( 2.39| ). The deviation from their exact counterparts 



( 2.11 ) and ( 2.12| ) is used as a measure for the goodness of the approximation. We pro- 



ceed by simply stating the results in Sects. 3.2, 3.3, 3.4 and leave the discussion of the 



consequences to Sect. 3.5 
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3.1 Four levels of coarse graining 

In this section, we introduce various ways to employ coarse graining, later used for each 



model in Sects. |3.2| , |3.3| and In all of the cases, the angular diameter distance is related 
to the luminosity distance by the trivial factor, d^iz) = (1 + z)'^dA{z), so we only give 
diiz) for the different cases. We list the cases here with their later referred names in bold: 

• No coarse graining — Exact LTB 

Here the observables are determined by the exact Eqs. ( 2.11| ) and ( |2.12| ), which we 



use as a baseline for the other cases. Conceptually, this model serves as a source of 
the actual observations a creature living in the exact LTB universe would make. 

A single averaging scale R — Single scale 

In this case, the observables are determined by the Eqs. ( p. 38 ) and ( p. 391) with R{z) = 
constant = R: 

4(z) = 2^0-1 (ii)(l + z)(l--^=L=) . (3.1) 



When considering a model that fits the supernova observations in Sect. we take 
the averaging scale R as the present-day physical distance to the object with the 
highest redshift in the supernova sample, R = rLTB(-Zmax)) numerically computed 
from Eqs. and (p.Sj). Other choices would at most correspond to different values 



of the Hubble constant, as is evident in Eq. (3T). For periodic inhomogeneities in 



Sect. |3.4| , we choose R = 27rro, where tq is the wavelength of the inhomogeneities. 

A different averaging scale R{z) for each redshift z — Running scale case. 

We further divide this into two distinct subcases, according to the explicit form of 
the function R{z): 

1. R{z) = f(z) — Running scale with averaged geodesies 

In this case, we take R{z) as the present-day physical distance to each redshift. 



determined by the averaged geodesies ( 2.37 ), in which we take R = ?'LTB('Zmax) or 



R = 2-KrQ. Choosing again R = f{z) would lead to the iterative use of Eq. (|2.37| ), 
perhaps ultimately converging to rLTB(-2) and making it no different from our 
next case. Although not used here, this could be a practical way of computing 
the distance in more realistic models where the exact result is unattainable. The 
observables are determined by 

4(.) = 2F„-' (2fl„-(ii)(l - ^)) (1 + - ^) . (3.2) 

R{z) = rLTB(^) — Running scale with exact geodesies 
Here we take the present-day physical distance to each object as the running 
averaging scale, determined by the geodesies of the exact metric, Eqs. ( p. 7]) and 
(|2.8D, yielding for the observables: 



dL{z) = 2H^^{nj,^{z)){l + z)[l-^^=^ . (3.3) 
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3.2 Analytic considerations — small z behavior 



There are two comphcations in comparing different realizations of the scale dependent 
averaging without specifying the boundary condition function HQ{r). Firstly, there does 
not exist expressions for the exact observables ( 2.1l| ) and ( |2.12| ) in terms of elementary 
functions. Secondly, already the relative comparison of the expressions for the averaged 
observables ( ^.Ij) , ( ^ ) and (|3.3| ) is unfeasible without knowing something about the func- 
tion Ho^r). 

A natural solution for the two problems is to calculate Taylor expansions for both 
the coarse grained observables of Eqs. ( |3.1| ), (|3.2[) , (3.3), and the exact expressions ( 2.11] ) 



and ( 2.12 ). It is straightforward to take the expansions to any desired order, but as the 
expressions become more complicated and unillustrative for higher orders, we give them to 
third order in redshift z. The price to pay is that the analytic comparison is valid only at 
small redshifts; we have numerically tested that the expansions are usable up to redshifts 
z ~ 0.2. 

We list the expansions of the luminosity distance here for the different cases following 



the entitling of Sect. 3.1: 



• Exact LTB 



8^3/^2(0)^^^4(0) 2Hi{0) 



1 



z' + 



1 ^ 1 if^(O) 



+ 2 



z^ + 0(z^ 



(3.4) 



R{z) = R — Single scale 



H-\R) 



Z+-Z'' 

4 



1 



z^ + 0{z^ 



(3.5) 



R{z) = r{z) — Running scale with averaged geodesies 



1 1 

+ 



^^(0) 



+ 



Ho{0)HoiR) 
^^(0) \ 



z' + 



8 2Ho{0)Ho{R) Hi{0)H^{R) 2 Ho{0)Hi{R) J 



z^ + Oiz'' 



(3.6) 



R{z) = rLTB(^) 



Running scale with exact geodesies 



i7^(0) 



z + 



1 1 H'JO) 



1 

4 ~ 
2Hl{0) 



z' + 



z'^ + Oiz^ 



(3.7) 
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3.3 Acceleration without backreaction 

The expansion of the dust dominated flat LTB universe can have neither local nor averaged 
acceleration but nevertheless can, as shown e.g. in Sect. 3.2 of Ref. [27], fit the supernova 
observations. Thus, from the observational point of view, the model can have effective 
acceleration. Since the backreaction vanishes in this model, the only possibility to account 
for the effect within the Buchert averaging formalism seems to be the running smoothing 
scale. 

In this section, we use the inhomogeneity profile found in Sect. 3.2 of Ref. [27], that 
gives a good fit to the Riess et. al. gold sample of 157 supernovae [1]. The boundary 
condition function of this model is given by 

Ho{r) = H + AHe-''/''o , (3.8) 

where the parameters have the values H + AH = 65.5 km/s/Mpc, AH = 16.8 km/s/Mpc 
and ro = 1400 Mpc. 

In Fig. |l], we plot the exact angular diameter distance of this model with the different 
averaged cases introduced in Sect. |3.1| . Instead of the dL{z) observed from the supernovae, 
we prefer to use dA{z) = dL{z)/{l + z)'^, as it is more slowly increasing function of z, 
making the differences between the various cases easier to detect. Moreover, the relative 
deviations {dA{z) — dAiz))/dA{z) are displayed in Figs. ^ and |3|, where dA{z) is the exact 
result of Eq. ( |2.12| ) and dA{z) stands for the averaged expressions ( p.lj ), (^), (|3.3|). Due to 
the the general relation di^z) = {1 + z)'^dA{z), the figures represent the relative deviations 
of the luminosity distance as well. 

Finally, we use the goodness of the fit to the Riess et. al. supernova data as an additional 
measure of the deviation, given by 

n=l ^ ' 

where an is the estimated error of the measured luminosity distance d'f^^{zn) to a source 
with redshift z„. In the different cases, diiz) of Eq. ( |3.9| ) is calculated from Eqs. ( 2.11| ), 
(|3.l| ), (^^), ( ^^ ) and gives the following values: 

• Exact LTB: 

= 1-12 (3.10) 

• Single scale 

with R = 4944 Mpc : = 4-35 (3.11) 
with R = 905 Mpc : = 2-06 (3.12) 

• Running scale with averaged geodesies: 

= 1-11 (3.13) 

• Running scale with exact geodesies: 

= 1-11 (3.14) 
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z 



O Exact LTB 


Single scale R - r^^^ ( 1 .755) 


■R(z)-T{z) — 


- R(z) ^r^^j^{z) 



Figure 1: The angular diameter distance dA{z) of the bubble model with Ho{r) = H + AHe 
where H + AH = 65.5 km/s/Mpc, AH = 16.8 km/s/Mpc and tq = 1400 Mpc. The black circles 
represent the exact LTB solution, whereas the blue, green and red curves correspond to the following 
coarse grained cases: Blue — a single averaging scale, chosen to be rLTB(1.755) — 4944 Mpc, where 
z = 1.755 is the redshift of the farthest supernova in the sample. Green — the averaged physical 
distance f{z) as the running smoothing scale. Red — the exact physical distance ruxBiz) as the 
running smoothing scale. 
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Single scale R - r^^^ ( 1 .755) — 


— R(z) =y(z) 


O R(z) =r^^^(z) 





Figure 2: The relative deviation {dA{z) — dA{z))/dA{z) = {dL{z) — d^iz)) / dL{z) of the averaged 
angular diameter (or luminosity) distance dA{z) from the exact value dA{z) for the bubble model 
of Sect. 3.3 in the following cases: Blue — a single averaging scale, chosen to be r ltb(1-755) = 
4944 Mpc, where z = 1.755 is the redshift of the farthest supernova in the sample. Green — the 
averaged physical distance f(z) as the running smoothing scale. Red — the exact physical distance 
fUT'B.iz) as the running smoothing scale. 
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— R(z) =y(z) 


O R(z) =r^^^(z) 







Figure 3: Same as in Figure 2, but for larger redshift range, z = 0...10. 



3.4 Periodic inhomogeneities as a toy model for structure 

Perhaps the closest representative of structure under the assumption of spherical symmetry 
is achieved with a periodic boundary condition function HQ(r) [33]. Hence, we take 

Ho{r) = H + AHsinr/ro , (3.15) 

where the values H = 65.5 km/s/Mpc, AH = 1.64 km/s/Mpc and tq = 95 Mpc have 
been chosen to make the plots as illustrative as possible. We do not consider more intense 
inhomogeneities in order to keep the relation tltbI-s) single- valued to the redshift z = 1. 

We plot the exact angular diameter distance of this model with the averaged results in 
Fig. 1^, and the relative deviation {dA{z) — dA{z)) / dA{z) in Fig. |5|, for each of the different 



cases introduced in Sect. 3.1 
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— Single scale R = 2 7i 


R(z)=7(z) — 


- R(z) ^r^^j^{z) 



Figure 4: The angular diameter distance dA{z) for the model of Sect. 3.4 with periodic inhomo- 
geneities, HQ{r) = H + AH s'm{r / tq) , where H = 65.5 km/s/Mpc, AH = 1.64 km/s/Mpc and 
To = 95 Mpc. The black circles represent the exact LTB solution, whereas the blue, green and red 
curves correspond to the following coarse grained cases: Blue — a single averaging scale, chosen 
to be one oscillation period 27rro = 598 Mpc. Green — the averaged physical distance f{z) as the 
running smoothing scale. Red — the exact physical distance ^ltbI^) as the running smoothing 
scale. 
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Single scale R = 2 7i r 



R(z) =y(z) 



R(z) 



LTB 



Figure 5: The relative deviation ((Ja(z) — dA{z))/dA{z) = [dL{z) — d^iz)) / dL{z) of the averaged 
angular diameter (or luminosity) distance dA{z) from the exact value dA{z) for the model of Sect. 
3.4 in the following cases: Blue — a single averaging scale, chosen to be one oscillation period 
27rro — 598 Mpc. Green — the averaged physical distance f{z) as the running smoothing scale. Red 
— the exact physical distance rijYB{z) as the running smoothing scale. 
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3.5 Discussion of the results 



Let us discuss here the results of Sects. 3^, 3^ and "iA by taking each coarse graining level 
of Sect. |3]l| into separate consideration: 



3.5.1 R{z) = R - Single scale 

In this case, already the general form of the luminosity distance ( |3.1[ ) reveals the essential 
point: averaging over a single scale is equivalent to using the perfectly homogeneous and 
flat matter dominated FRW model. The freedom to choose the averaging scale R only 
corresponds to fixing the value of the effective Hubble constant Ho{R) in Eq. (p.ll). 

The supernova data fit in the bubble model of Sect. shows that when averaging over 
the domain that contains all the supernovae in the sample {R = 4944 Mpc), the resulting 
= 4.35 represents a huge deviation from the exact result, = 1.12. We get a better fit 
by choosing R = 905 Mpc, because this gives the effective Hubble constant in Eq. its 
best fit value for the fiat matter dominated FRW case, i^o(905 Mpc) = 57.5 km/s/Mpc. 
However, even this value gives = 2.06, which is still a way too large deviation from the 
actual result. Besides, we have no a priori physical justification to pick up the particular 
averaging scale R = 905 Mpc. Overall, it is clear that no single averaging scale can give 
an acceptable approximation for the bubble model. 

Perhaps the most interesting feature of the single scale case becomes evident in Figs. 
Q and ^: even though the model of Sect. 3.4 with periodic inhomogeneities is homogeneous 



on large scales, the single averaging scale still leads to unwanted deviations. Whether it is 
an artifact of the employed spherical symmetry, with light inevitably propagating through 
all the layers of structure, or a more general phenomenon, remains an open question. 

Altogether, the expression of the averaged luminosity distance ( |3.lD , along with its 
Taylor expansion (|3.5D , and the figures || - ^ all confirm the conclusion that averaging 
over a single scale gives a too coarse-grained description at least for the fiat LTB universe. 
The inadequacy of the averaging procedure to account for the observations in the LTB 
universe was already suggested in Ref. [27], but our results bring out the essential point 
that the conclusion is valid only under the assumption of a single smoothing scale, as we 
next discuss. 

3.5.2 R{z) = f{z) — Running scale vi^ith averaged geodesies 

It has been speculated that averaging cosmological inhomogeneities would be useful only 
when applied to systems with statistically homogeneous distribution of small scale irregu- 
larities [19]. Nevertheless, our results indicate that it is possible to exploit averaging also 
for large scale inhomogeneities, if the single smoothing scale R is promoted to a redshift 
dependent function R{z). 

On physical grounds, one could have expected some improvement in accuracy of the 
observables when the single averaging scale R is replaced by the averaged present-day 
physical distance ( |2.37| ) to each object at redshift z. However, the amount of precision 
achieved with this generalization is both surprising and a very welcome result. Indeed, 
the congruence between this approximation and the exact results is evident in all of the 



comparisons made in Sects. 3.2, 3.S and 3.4 as we next specify in more detail. 



-17- 



When comparing the Taylor expanded luminosity distance (|3.q ) with the expansion of 
the exact LTB case (3.4), one sees that to second order the results are almost identical. 
They become exactly identical if we take -R = as the averaging scale for the geodesies, 
corresponding to the use of the local Hubble parameter Hq{0) in Eq. ( p. 37 ), or if Hq{R) = 
Hq{0), as is the case for periodic inhomogeneities of Sect. |3.4| . Moreover, even the third 
order terms carry the same functional dependence on HQ{r) as the exact case of Eq. (p.4|), 
albeit with slightly different prefactors. 

The supernova data fit in the bubble model of Sect. 3^ illustrates the power of the 
running scale approach in practical applications. When employing the running scale, the 
goodness of the fit changes from the puny = 4.35 of the single scale case to the excellent 
fit = 1.11, which is within one percent of the correct result = 1.12; the improvement 
is manifest in Figs. |l] and |2| as well. The result also elucidates how the scale dependence 
can produce apparent acceleration even in the absence of backreaction; for a thorough 
discussion of the apparent acceleration, see page 9 of Ref. [27]. Anyhow, in the real 
universe, one could expect the effective acceleration to arise from the interplay between 
scale dependence and backreaction, but the possibility remains that either — or neither for 
that matter — of them plays the dominant role. 

The toy model of structure in Sect. |3.4| also manifests the advantage of the running 
smoothing scale: it is evident in Figs. ^ and ^ that the running scale follows the oscillations 
of the exact observables whereas the single scale case simply fails to. 

Problems arise, when going beyond the supernova fits to higher redshifts. Indeed, Fig. 
m reveals that at z > 2 the running scale falls short of the 0(1%) accuracy compared to 
the exact observables. There is a plausible physical explanation for this: in the employed 
LTB model, the growth of inhomogeneities backwards in time^ makes them more important 
at higher redshifts and cannot be encapsulated in the present-day spatial averages. One 
can still argue that in a more realistic model the problem would be alleviated, since the 
inhomogeneities of the real universe are expected to grow forwards in time. Overall, we 
suppose the averaging with the running scale would be conceivable at least up to z ~ 2; for 
higher redshifts, one could then use coarser approximations, such as the perturbed FRW 
models, presuming inhomogeneities really were of less importance in the past. 

It would, nevertheless, be desirable to remedy the problems at higher redshifts. For 
this, we have considered various forms for the running smoothing scale R{z), beyond the 
ones introduced in Sect. ^T^. The outcome is that there seems to be no functions R{z) that 
would both have a physical basis, and improve the approximation over the whole range of 
redshifts. In particular, it seems difficult to allow for the time evolution of the universe 
within the running averaging scale approach. Perhaps a better way to take into account the 
time evolution would be to average over the past light cone. However, this would require 
a complete revision of the basic formalism. 

Finally, a great virtue of using the geodesies of the averaged metric is that one does 
not need to solve the exact geodesic equations. Moreover, one can improve the accuracy 
of this method by iterative use of the averaged distance-redshift relation (2.37), as argued 
in the paragraph before Eq. (|3.2D. 



^See Eqs. (2.5), (|2.6[) and the paragraph thereafter. 
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3.5.3 R{z) = rijTB{z) — Running scale with exact geodesies 

There are two conceptual steps in coarse graining needed to calculate the observables 
[32]. Firstly, the step from the exact Einstein equations to the Buchert equations, and 



secondly, the step from the exact metric (2.1) to the average metric (2.34). To quantify the 



approximation of the latter step, we studied a case where only the field equations have been 
averaged, but instead of the average metric, the exact metric determines the geodesies. 

The outcome is that there are only minor deviations in the observables between the use 
of the averaged and the exact geodesies. This becomes apparent in the congruence between 
the Taylor expansions of the luminosity distance ( |3.6| ) and ( p.?! ), between the resulting 
for the supernova data fit in Eqs. ( |3.13| ) and (3.14), and between the red and green curves 



in Figs. 1^ - ||- Although corresponding only to a slight correction, it is still evident in all of 
these comparisons that the exact physical distance gives the most accurate approximation. 
Anyway, due to the good congruence of the results between the averaged and the exact 
geodesies, the feasibility becomes the deciding factor. Indeed, in more realistic models 
of the universe, the exact geodesies are beyond computation. Overall, perhaps the best 
solution in practice is to use the averaged geodesies and, if needed, use Eq. (|23^ ) iteratively 



as explained in Sect. 
4. Conclusions 

We have considered the role of scale dependence in the Buchert averaging formalism, using 
the spherically symmetric and spatially flat LTB dust universe as a testing ground. The 
vanishing of both the backreaction and the spatial curvature scalar ^^^R makes this an 
ideal model to capture the effect of the averaging scale on the observables, because then 
the only difference from the flat FRW model is the explicit dependence on the smoothing 
scale. From the resulting Buchert equations ( p. 31 ), ( p.32|) , ( p.33|) and the average metric 



(2.34), we have derived the luminosity and angular diameter distances ( p. 38 ) and ( 2.3S| ) 



carrying the scale dependence as well. By employing a redshift dependent averaging scale 
R{z), we have compared these observables to the exact expressions (2.11) and ( p. 12 ). The 



physical reason to use a different averaging scale for each object at redshift z is clear: the 
distance the observed light propagates depends on how far the object is. 

Our principal result is that when the conventional single averaging scale is replaced by 
the physical distance to each object at redshift z, the relations ( 2.38| ) and ( p. 39]) become 



significantly closer to the actual observables without complicated computations. Although 
some improvement could be expected on physical grounds, the 0(1%) precision at z < 
2 makes the result a welcome surprise. Indeed, contrary to some previous speculations 
[19,27], the result suggests that the averaging procedure can be exploited for large scale 
inhomogeneities, presuming the running smoothing scale is employed. Although considered 
merely under the assumption of spherical symmetry, we expect the running scale to show 
its full advantage only when applied to more irregular large scale inhomogeneities, such as 
the recently observed voids [21, 22]. Naturally, the running scale can also be applied even if 
a single scale would suffice, since then it reduces to give the same predictions as the single 
scale approach. 
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By Taylor expanding the averaged observables and their exact counterparts, we have 
demonstrated the increase in accuracy for generic inhomogeneities up to redshifts z ~ 0.2. 
In addition, we have numerically confirmed it for redshifts up to ^ ~ 2 using two explicit 
inhomogeneity profiles: a bubble inhomogeneity that fits the supernova observations and 
periodic inhomogeneities as a toy model for structure. Within the bubble model, we found 
that the running smoothing scale can account for the apparent acceleration even with- 
out backreaction. Consequently, it could be at least as important as the backreaction in 
explaining dark energy as an inhomogeneity induced illusion. 

At redshifts z > 2, the method breaks down. The plausible physical explanation for 
this is that the present-day spatial averages cannot capture the time evolution of the exact 
model. Whether the problem can be solved within the already developed formalism, or a 
generalized approach such as averaging over the light cone would be needed, is something 
we hope to address in a future work. In any case, the time evolution of the fiat LTB model is 
contrary to the observed structure formation, since the inhomogeneities of the model grow 
towards large redshifts. Hence, perhaps in the real universe one could apply the running 
scale in averaging at low redshifts, and at higher redshifts, resort to the conventional 
perturbed FRW models. 

We have employed two ways to compute the present-day physical distance R{z) to ob- 
jects with redshift z: firstly, using the geodesies of the average metric and, for comparison, 
using the geodesies of the exact metric. The outcome is that using the exact geodesies 
gives only a marginal correction compared to the averaged geodesies and, requiring the 
exact solution, would in general lead to complicated computations. Overall, maybe the 
best method in practice would be an iterative use of the averaged distance-redshift relation 



( 2.37 ), as explained in Sect. |3.1| . 

Finally, there are problems in cosmological averaging we have not addressed in this 
work (see [19,23,25,34-38]). Perhaps the most notable one within the Buchert formal- 
ism is that the averaged equations ( |2.13 ), ( 2.14| ) and (|2.15| ) contain three equations for 
four unknowns. Therefore, more information is in general needed to solve the equations. 
Whether this makes the whole approach impractical in calculating observables outside the 
exact solutions of the Einstein equations, is an open question. 
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